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Abstract 

A  numerical  simulation  based  on  the  optimal  control  theory  was  developed  to  assess  the  performances,  dimensions  and  weight 
characteristics  of  skutterudite-based  solar  thermoelectric  generators  (STG)  for  satellite  missions  at  distances  close  to  the  Sun.  It  is  shown 
that  the  STG  designs  based  on  these  advanced  thermoelectric  materials  are  potentially  attractive  for  this  type  of  application  and  could 
advantageously  replace  the  standard  power  sources  such  as  the  solar  (cells)  panel  or  the  radioisotope  thermoelectric  generators. 

©  2002  Elsevier  Science  B.V.  All  rights  reserved. 

Keywords:  Thermoelectric  devices;  Thermoelectric  energy  conversion;  Skutterudites;  Optimal  control  theory 


1.  Introduction 

So  far,  the  primary  power  sources  used  on  unmanned 
spacecraft  applications  are  the  solar  arrays  and  the  radio¬ 
isotope  thermoelectric  generators  (RTG).  However,  when 
working  at  distances  close  to  the  Sun,  the  use  of  solar  cells  is 
limited  by  two  factors:  (1)  the  high  incident  solar  heat  flux 
tends  to  produce  high  solar  cell  temperatures  and  (2)  the  high 
incident  radiation  flux  causes  cell  damage.  Both  factors 
yield  rapid  deterioration  of  the  solar  cell  performance.  Even 
if  the  solar  cells  are  temperature  controlled  and  protected 
against  the  radiation  damages,  they  become  less  attractive  in 
near-Sun  missions  because  of  increased  weight  and 
decreased  reliability.  RTG  have  been  extensively  used  on 
interplanetary  US  satellites,  but  are  not  used  since  more  than 
10  years  on  board  European  scientific  satellites.  The  envir¬ 
onmental  risks  in  case  of  explosion  during  the  ascent  and  the 
difficulty  to  import  such  system  were  the  main  reasons  to 
abandon  this  technology. 

Among  candidates  other  than  RTG  that  could  operate 
efficiently  and  power  satellites  flying  near  to  the  Sun,  the 
solar  thermoelectric  generators  (STG)  could  be  an  interest¬ 
ing  alternative  solution,  because  the  thermoelectric  materi¬ 
als  are  in  general  less  sensitive  to  both  the  temperature  and 
the  radiation  damages.  Systems  based  on  STG  have  already 
been  described  in  the  past  [1-7]. 


*  Corresponding  author.  Tel.:  +33-3-83-58-41-63; 
fax:  +33-3-83-57-97-94. 

E-mail  address:  bertrand.lenoir@mines.inpl-nancy.fr  (B.  Lenoir). 


A  complete  STG  consists  of  the  collector/concentrator, 
the  thermopiles,  the  insulation  and  the  radiator.  The  material 
selection  for  each  of  these  four  parts  is  of  the  greatest 
importance  in  the  overall  performance  and  weight  charac¬ 
teristics  of  the  system  and  shall  be  considered  in  details. 
Referring  to  the  approximate  relation: 

r]~\ZAT  (1) 

it  is  easily  seen  that  in  order  to  achieve  a  high  efficiency  r\, 
the  temperature  difference  AT  across  the  thermopile  should 
be  as  high  as  possible,  and  the  thermocouples  constituting 
the  thermopile  should  possess  a  high  figure-of-merit  Z. 
Because  Z  of  established  thermoelectric  materials  differ, 
as  do  the  ranges  of  temperature  difference  across  which 
they  can  be  operated,  it  is  normally  found  that  most 
thermoelectric  materials  yield  performance  values  that  differ 
somewhat  from  each  other.  Over  the  past  few  years,  much 
effort  has  been  devoted  to  achieve  high  Z  values  through 
the  exploration  of  new  classes  of  materials  such  as  semi¬ 
conductors  and  semimetals  with  complex  open  crystal 
structures  allowing  for  very  low  thermal  conductivity 
[8,9].  Among  them,  skutterudites  based  on  CoSb3  have 
been  recently  promoted  as  possible  new  thermoelectric 
materials  [8-10].  n-  and  p-type  materials  have  been  synthe¬ 
sized  and  were  found  to  have  higher  figure-of-merits 
than  those  of  the  established  Si-Ge  based  alloys  in  the 
temperature  range  300-1000  K. 

In  order  to  assess  the  performance,  dimensions  and  weight 
characteristics  of  a  skutterudite-based  STG  for  near-Sun 
missions,  a  physical  model  has  to  be  developed.  The  number 
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of  variables  is  important  and  it  is  necessary  to  solve  a 
number  of  optimizing  problems.  The  algorithm  for  the 
optimization  was  generally  based  on  classical  energy  bal¬ 
ance  equations,  which  did  not  take  into  account  neither  the 
dependences  of  the  thermoelectric  parameters  on  the  tem¬ 
perature,  nor  the  contact  resistances  which  may  result  in 
great  errors  [1,6].  The  above-mentioned  approximations  can 
be  avoided  using  modern  methods  of  optimization,  for 
instance  the  optimal  control  theory  [11,12]. 

In  this  paper,  we  present  a  numerical  simulation  using  the 
optimal  control  theory  to  design  a  STG  based  on  skutterudite 
thermoelectric  materials.  In  a  first  place,  the  mathematical 
development  is  presented.  It  is  then  used  to  study  parame¬ 
trically  the  performance  and  the  mass  analysis  of  three 
different  STG  concepts.  These  results  are  discussed  in 
details  in  relation  with  the  STG  design.  The  imposed 
requirements  for  the  STG  were:  (i)  the  orientation  perpen¬ 
dicular  to  the  Sun,  (ii)  the  production  of  400  W  electrical 
power  and  (iii)  the  working  distance  of  0.45  a.u.  from  the 
Sun  that  corresponds  to  Mercury  (average). 


2.  Design  method  of  STG 

2.7.  Physical  model  for  STG 

A  complete  STG  consists  of  a  solar  concentrator  or 
collector,  a  set  of  thermopile  units,  insulating  parts  and  a 
radiator  (Fig.  la).  Several  types  of  concentrators  can  be  used 
for  increasing  the  solar  intensity  and  acquiring  an  increase  in 
thermal  flux  on  the  hot  junction  surface  of  the  thermopiles. 
The  configurations  of  interest  discussed  in  earlier  papers 
involve  a  flat-plate  collector  [1-3]  or  a  compact  design  with 
a  single  concentrator  or  with  multiple  solar  mini-concen¬ 
trators  [4,5].  A  single  thermopile  consists  of  several  thermo¬ 
electric  couples  (n-  and  p-type),  connected  usually  in  series, 
sandwiched  between  two  ceramic  plates  and  sealed  around 
the  edges  (Fig.  lb).  This  structure  allows  to  connect  the 
thermopiles  in  parallel  within  the  STG  providing  more 
reliable  operation  of  the  thermoelectric  device  as  a  whole. 
The  space  between  the  collector  and  radiator  plates  not 
occupied  by  the  thermopiles  is  filled  with  a  thermal  insulator 
to  reduce  the  shunt  heat  losses. 

The  schematic  diagram  of  a  flat-plate  STG,  shown  in  Fig.  1, 
is  used  to  establish  the  appropriate  mathematical  relationships 
of  the  STG  physical  model.  For  simplicity,  the  heat  flow  is 
assumed  to  be  one-dimensional.  This  assumption  is  reason¬ 
able,  since  the  collector  and  radiator  plates  are  usually  thin 
and  made  of  materials  with  good  thermal  conductivity.  There¬ 
fore,  the  heat  flow  along  the  STG  due  to  local  temperature 
gradients  is  neglected  as  it  is  usually  compared  to  the  heat  flux 
across  the  STG.  The  heat  losses  from  the  edges  of  the 
insulator,  collector  and  radiator  plates  are  neglected. 

The  collector  plate  is  assumed  to  be  exposed  to  a  per¬ 
pendicular  solar  flux.  The  heat  absorbed,  Qa ,  is  basically 
dependent  on  the  distance  from  the  Sun  d  and  can  be 


insulator 

radiator 


Fig.  1.  (a)  The  principle  scheme  of  a  flat-plate  STG.  (b)  Schematic 
diagram  of  a  single  thermoelectric  couple  constituted  of  n-and  p-type 
semiconducting  branches. 


determined  as  follows: 


Ga  = 


$a  Sc 


Al 


(2) 


where  aa  is  the  absorptivity  of  the  collector  coating,  Sc  the 
solar  constant  at  1  a.u.  in  space  (Sc  =  1370  W/m2)  and  Als 
the  collector  area.  The  amount  of  heat  Qe  re-radiated  from 
the  collector  surface  into  space  at  0  K  (assumed  tempera¬ 
ture)  is  given  by: 

2e  =  SiVSbAiTi  (3) 

where  Si  is  the  emissivity  of  the  collector  coating,  crSB  the 
Stefan-Boltzmann  constant  and  Tx  the  temperature  of  the 
collector  surface. 

The  shadow  face,  i.e.  the  radiator,  is  cooled  by  the 
radiation  of  heat  flowed  through  the  thermopiles  and  the 
insulation  to  deep  space.  This  amount  of  heat,  Qr,  lost  into 
space  is  given  by: 


<2r  =  (4) 

where  s2  is  the  emissivity  of  the  radiator  coating  and  A2,  T4  the 
area  and  the  temperature  surface  of  the  radiator,  respectively. 

In  operation,  the  total  heat  transported  through  the  STG  is 
the  amount  of  heat  absorbed  by  the  collector  from  the  Sun 
minus  the  energy  converted  to  electric  power. 

Written  in  terms  of  heat  flux  densities  for  steady-state 
conditions,  the  previous  energy  balance  equations  become: 


q a  qe  -  „  .  (Ti  t2) 

Kitli 

(5) 

1  , 

qr  =  — —  (73  —  74) 

H  R2h2  } 

(6) 

H.  Scherrer  et  al./  Journal  of  Power  Sources  115  (2003)  141-148 


143 


q,K}  -  qtK\  -  qh{T2 ,  T3)  -  qfK2  -  1)  =  0  (7) 

qvK2  -  q i [K2  -  1)  -  qc(T2,  T3)  =  0  (8) 


where  R\  and  R2  are  the  specific  heat  resistances  of  collector 
and  radiator  materials,  h\,  h2  the  thicknesses  of  these 
materials,  K\  =  (Ai/Ate),  K2  =  (A2/ATE)  the  ratios  of 
collector  Ax  and  radiator  A2  areas  to  the  thermopiles  area 
Ate. 

In  the  set  of  Eq.  (5)-(8),  q a  represents  the  density  of  heat 
absorbed  by  the  collector  from  the  Sun: 


aa5c 


(9) 


q&  denotes  the  density  of  heat  re-radiated  by  the  collector 
surface: 


<7e  =  ~  =  ei  r7SBrf  (10) 

qr ,  is  the  density  of  heat  emitted  by  the  radiator  surface: 

qr  =  ^  =  e2^SB?4  (11) 

^2 

and  qi  is  the  density  of  heat  transported  through  the  insula¬ 
tion: 


qi  = 


Qi 

a2  —  ate 


A 

h 


(T2  -  t3) 


(12) 


where  X  is  the  thermal  conductivity  of  the  insulator  and  the 
thickness  of  insulation  between  the  collector  and  the  radia¬ 
tor. 

The  densities  of  heat  absorbed  at  the  hot  junctions  of  the 
thermocouples  and  rejected  at  the  cold  thermoelectric  junc¬ 
tions  are  denoted  as  qh,  qc ,  respectively.  These  two  heat  flux 
densities  depend  on  the  temperatures  of  hot  T2  and  cold  T3 
sides  of  the  thermopiles  (see  Fig.  1).  They  are  also  assigned 
by  the  operating  mode  of  the  thermocouples.  Expressions 
(5)— (8)  with  the  additional  relationships  for  calculation  of 
heat  densities  qh  and  qc  form  a  system  of  four  nonlinear 
equations  which  for  any  couple  of  given  K\  and  K2  values 
define  uniquely  the  unknown  temperatures  7\,  T2,  T3,  T4. 
Using  these  data,  all  unknown  characteristics  of  the  STG  can 
be  estimated. 

The  relationships  for  heat  densities  qh  and  qc  can  be 
obtained  from  the  theory  of  thermoelectric  conversion,  using 
the  detailed  optimization  of  thermoelectric  device  perfor¬ 
mance. 


2.2.  Design  optimization  of  thermoelectric  devices 

Several  methods  have  been  used  for  designing  and  opti¬ 
mizing  thermoelectric  power  generators  with  solar  heat 
source  [1,2, 6, 7].  The  main  imperfection  arises  from  the  fact 
that  most  of  them  are  based  on  the  relationships  of  heat 
balance  at  the  hot  and  the  cold  junctions  of  the  thermoele¬ 
ments  which  are  only  valid  when  Joule  and  Thomson  heats 
are  small  in  comparison  to  the  conduction  heat  and  when 


electrical  resistivity,  thermal  conductivity,  and  Thomson 
coefficient  are  not  temperature  dependent.  In  most  of  the 
earlier  published  papers,  Thomson  heat  is  even  neglected 
[1,6].  This  approach  is  appropriate  only  for  materials  with 
little  variation  of  the  Seebeck  coefficient  with  temperature. 
But  in  the  case  of  thermoelectric  generation  under  large 
temperature  gradients,  the  dependence  of  the  material  para¬ 
meters  on  temperature  is  appreciable  and  should  be  taken 
into  account.  Moreover,  some  important  factors  such  as 
thermal  and  electrical  contact  resistances  are  not  always 
considered.  It  is  often  assumed  that  thermocouples  have  a 
fixed  configuration  and  optimization  is  made  by  calculation 
of  the  device  performance  for  combination  of  values  for 
some  variables  [2,3].  In  such  methods  the  number  of  com¬ 
binations  is  limited,  which  may  lead  to  inaccurate  optimiza¬ 
tion  results. 

More  recently,  a  novel  numerical  method  suitable  for 
computer  implementation  has  been  developed  for  design 
and  optimization  of  performance  of  thermoelectric  systems. 
It  enables  to  overcome  the  shortcomings  of  the  previous 
methods.  This  numerical  simulation,  based  on  the  mathe¬ 
matical  theory  of  optimal  control  [11,12]  is  very  attractive 
for  the  optimization  and  the  prediction  of  the  thermoelectric 
module  performance.  It  has  been  developed  in  details  for 
thermoelectric  coolers  and  has  been  further  extended  suc¬ 
cessfully  to  optimize  thermoelectric  generators  [13,14].  We 
propose  here  to  apply  this  method  to  extract  the  heat  flux 
densities  qh  and  qc  at  the  hot  and  cold  junctions  of  the 
thermopiles  (T2  and  T3)  involved  in  the  system  of  balance 
Eqs.  (5)-(8). 

The  initial  data  for  the  STG  design  are  the  output  elec¬ 
trical  power  WL  and  voltage  VL  on  the  external  load.  It  is  also 
important  to  select  the  mode  of  operation  of  the  generator. 
For  this  purpose,  the  STG  efficiency  can  be  written  as 
follows: 


Wl  w 

,STa  =  a=i^ 


(13) 


where  w  =  (  Wl/Ate)  is  the  output  electric  power  per  unit  of 
thermoelement  area. 

In  our  case  the  solar  heat  density  qa  is  fixed.  Therefore,  for 
any  given  K\  value,  the  maximum  STG  efficiency  mode  is 
equivalent  to  the  maximum  output  power  density  w,  which 
is: 


w  =  qh~qc 


(14) 


The  problem  is  then  to  maximize  the  functional  w. 

Temperature  T{x)  and  heat  flux  q(x)  distributions  within 
one  of  the  legs  (n  or  p)  of  a  thermocouple  at  a  distance  v  from 
the  cold  junction  (Fig.  1)  are  determined  by  the  system  of 
four  differential  equations  of  the  form  [13]: 


dT 

dx 

dq 

dx 


_^lT  _q 


K 


K 


2-2  •  -2 

U  J  T  ,  a/  ,r 

- T  H - q  H — 

K  K  ( 7 


> 

/  n,p 


(15) 
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where  a,  a,  k  are  the  Seebeck  coefficient,  the  electrical 
conductivity  and  the  thermal  conductivity,  respectively  and  j 
the  current  density. 

In  this  system,  the  thermoelectric  material  properties  a(7), 
(7(7),  k(T)  are  temperature  dependent  and,  approximated  by 
polynomials.  The  solution  of  Eq.  (15)  under  the  following 
boundary  conditions: 

Tn{  o)  =  7p(0)  =  r3,  Tn(l)  =  7p(/)  =  T2  (16) 

gives  the  opportunity  to  calculate  the  heat  flux  densities  at 
the  hot  and  cold  thermoelements  junctions  qh  and  qc  using 
the  expressions: 


For  such  a  calculation,  an  initial  value  must  be  given  to  all 
materials  properties  as  well  as  the  length  of  the  thermoele¬ 
ment  legs.  The  result  of  the  calculation  gives  a  complete 
spectrum  of  temperatures  T1?  . . .,  T4  and  correspondent 
values  of  maximum  output  power  density  w,  optimal  heat 
fluxes  gh,  qc ,  and  current  densities  jn,  jp.  Then,  for  the 
required  power  and  voltage  VL  on  the  external  load, 
the  complete  performance  and  parameters  of  the  STG 
structure  can  be  estimated  using  the  following  relationships: 

The  thermoelectric  generator  efficiency  q: 


qh  qh 


4h  =  +j2r c] 

If  (17) 

=  -yy<?(°)  -fr c] 

n,p 

where  rc  is  the  contact  resistance  and  /  the  thermoelement 
legs  height. 

The  relation  (17)  shows  that  qc  and  qh  depend  on  the 
current  densities  jn  p.  These  parameters  have  to  be  evaluated 
in  order  to  maximize  the  functional  w. 

The  optimal  control  method  allows  to  make  a  purposeful 
search  for  jn  p  values  [11,13].  According  to  this  method, 
optimal  current  densities  should  satisfy  the  following  rela¬ 
tionships: 


dw  |  rl  dHji^qJnJp)  dv  () 
9ja,p  Jo  <%,p 


(18) 


The  thermopiles  area: 


(21) 


w  qh  -  qc 

The  collector  surface  area  Ax  and  the  radiator  surface  area 

A2- 


Ai  =  TTiAte,  A2  =  K2Aje  (22) 

The  heat  fluxes  through  the  collector  plate  Qs  and  radiator 
plate  Qr: 

Qs  =  A\{cjdi  q&)i  Qr  =  A2qYl  (23) 

The  single  thermopile  area  (for  a  given  number  of  thermo¬ 
piles  AO: 


A 


Ate 


(24) 


where  the  Hamilton  function  H  has  the  form: 

h = + ^2)  (i9) 

n,p 

where  (/j,  /2)n, p  are  the  right  parts  of  (15)  and 
i/j  =  i/j2) np  is  the  pulse  vector  conjugated  to  the  phase 

variable  vector  y  =  (7)  q)np.  Its  determination  is  detailed 
elsewhere  [13]. 

The  numerical  method  of  successive  approximations  can 
be  recommended  to  solve  this  optimal  control  problem.  It 
allows  to  calculate  the  optimal  current  densities  jUrP  and  heat 
fluxes  qh ,  qc  values  respective  to  the  maximum  w. 

2.3.  Computer  algorithm  for  the  STG  design 


The  single  thermopile  electric  power  W,  voltage  V  and 
current  7: 


WL 

w  =  — - 

N 

(25) 

< 

II 

(for  series  thermopile  connection) , 

(26) 

II 

(for  parallel  connection) , 

(27) 

II 

(28) 

The  cross-sections  of  n-  and  p-type  legs  An?p: 


A„,p  =  2-  (29) 

7n,p 


The  algorithm  of  the  STG  design  method  described  above 
is  based  on  the  numerical  solution  of  the  heat  balance  system 
(5)— (8)  relative  to  the  unknown  temperatures  7),  . . .,  T4  for 
given  factors  K\  and  K2.  The  method  of  Newton’s  iterations 
has  been  used  for  solving  the  nonlinear  equations  system.  A 
subroutine  was  used  to  calculate  the  optimal  heat  densities 
qh  and  qc.  It  should  be  noted  that  the  equation  system  (5)— (8) 
includes  only  two  independent  nonlinear  Eqs.  (7)  and  (8) 
with  respect  to  the  independent  variables  7),  T4.  The  tem¬ 
peratures  T2,  T3  are  defined  in  Eqs.  (5)  and  (6). 


The  number  of  thermocouples  connected  in  series  per 
thermopile  ATE: 


Ate 


A 

(An  +  Ap) 


(30) 


The  weight  of  the  STG  can  be  directly  calculated  from  the 
dimensions  and  the  material  characteristics  used  for  the 
collector,  the  radiator,  the  insulation  and  the  thermopiles. 

As  a  result,  the  optimal  STG  performance  in  maximum 
conversion  efficiency  mode  and  simultaneously  in  maximum 
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power  output  per  unit  surface  of  thermoelectric  device  mode 
are  known  for  any  given  factor  K\  and  K2.  The  optimization  of 
the  performance  with  respect  to  Kx  and  K2  should  be  made 
through  the  variation  of  these  parameters. 

Based  on  this  algorithm,  a  computer  program  for  STG 
design  has  been  developed. 

2.4.  STG  concepts  and  material  selection 

The  proposed  mathematical  model  has  been  successfully 
used  for  the  design  of  three  STG  concepts  with  regard  to 
factors  Kx  and  K2.  In  the  first  concept,  collector  and  radiator 
plates  are  assumed  to  have  identical  areas  ( K\  =  K2).  It  is  a 
so-called  flat-plate  STG  (Fig.  1).  The  second  one  is  a  STG 
having  the  radiator  plate  area  larger  than  the  collector  plate 
area  ( K2  >  K\).  The  third  concept  is  concerned  with  a  STG 
configuration  with  a  single  solar  concentrator  or  with  mini¬ 
solar  concentrators  for  each  thermopile  unit,  in  place  of  the 
flat-plate  collector.  This  configuration  is  exemplified  in 
Fig.  2.  In  this  case  Kx  is  interpreted  as  a  concentration  ratio 
(usually  K\  >  K2).  The  advantages  of  one  STG  concept  on 
another  can  only  be  estimated  on  the  basis  of  the  calculation 
results. 

The  evaluation  of  the  performance  and  weight  character¬ 
istics  of  the  three  STG  designs  has  been  accomplished  for  a 
0.45  a.u.  operation  distance,  a  400  W  load  and  according  to 
the  following  material  selection.  As  previously  stated,  we 
propose  to  use  skutterudite  materials  as  thermoelectrics. 
This  novel  family  of  materials  has  received  considerable 
attention  this  last  decade  and  continues  to  be  extensively 
investigated  due  to  their  unique  electrical  and  thermal 
properties.  Their  use  in  thermoelectric  devices  for  terrestrial 
or  space  power  generations  has  also  been  considered  and  the 
results  seem  to  be  very  promising  [15].  Among  skutterudite 
materials,  filled  CoSb3  (p- type)  and  doped  CoSb3  (n-type) 
are  known  to  possess  high  thermoelectric  performance  in  the 
temperature  range  400-1000  K  making  them  very  attractive 
for  STG  applications.  The  density  of  these  materials,  rela¬ 
tively  low,  is  near  7.65  g/cm3. 

The  choice  of  the  collector  and  radiator  materials  was 
principally  guided  by  two  rules:  a  high  thermal  conductivity 
and  a  low  weight,  compatible  with  the  spacecraft  technol¬ 
ogy.  For  the  calculations,  we  tentatively  selected  collector 
and  radiator  of  aluminum  nitride  (AIN)  that  has  a  thermal 
conductivity  of  about  100  W/mK  at  high  temperature  and  a 


Fig.  2.  Schematic  diagram  of  a  STG  with  multiple  solar  concentrators. 


density  of  about  3.3  g/cm3.  It  is  assumed  that  the  collector  is 
covered  with  a  particular  coating  ensuring  high  solar  absorp- 
tance  and  low  emittance  (collector  absorptivity  aa  =  0.9; 
collector  emissivity  =  0.08).  Conversely,  the  radiator  is 
covered  with  an  optical  coating  ensuring  high  heat  rejection 
(radiator  emissivity  s2  =  0.75). 

The  material  ensuring  the  thermal  insulation  between  the 
collector  and  the  radiator  plates  must  have  a  thermal  con¬ 
ductivity  as  low  as  possible  in  order  to  reduce  shunt  heat 
transfer.  Our  calculations  were  performed  assuming  the  use 
of  a  thermal  insulator  from  Microtherm  Society,  having  a 
thermal  conductivity  of  0.035  W/mK  within  the  operating 
temperature  range  and  a  density  of  0.24  g/cm3.  This  insu¬ 
lator  is  located  between  the  individual  thermopiles  within 
the  STG,  so  that  the  minimum  height  of  insulation  to 
consider  is  the  height  of  the  thermopile.  Vacuum  is  assumed 
to  be  in-between  the  legs  of  the  thermopile.  The  value  of  the 
contact  resistance  between  the  cold  and  hot  electrodes  and 
the  legs  of  the  thermoelements  was  assumed  to  be  equal  to 
5  x  10-6  Q  cm2. 


3.  Results  and  discussion 

The  next  paragraphs  present  the  results  obtained  from  the 
numerical  simulation  for  the  three  STG  concepts.  The 
analysis  focuses  on  the  influence  of  Kx  and  K2  values  on 
the  dimensions,  performances  and  weight  characteristics  of 
the  STG.  The  complete  description  of  the  system  will  not  be 
described  here,  as  for  instance  the  number  of  thermocouples 
per  thermopile,  the  cross-section  of  the  n-  and  p-legs  and  the 
current  used.  These  data  can  be  estimated  in  a  straightfor¬ 
ward  way  once  N  and  Fare  fixed  using  the  set  of  Eqs.  (24)- 
(30)  exposed  in  Section  1.3. 

3.1.  Concept  1:  collector  and  radiator  of  same  area 

The  first  calculations  have  been  carried  out  for  a  flat-plate 
STG  model  with  collector  and  radiator  of  equal  area  A.  The 
modeling  results  obtained  for  thermoelectric  legs  of  1  cm  in 
height,  an  insulation  thickness  of  1.5  cm  and  plates  of 
collector  and  radiator  of  1  mm  in  thickness  are  summarized 
in  Table  1.  The  temperature  drops  AT  between  the  collector 
and  the  radiator  surfaces  (7\-r4),  as  well  as  between  the  hot 
side  and  the  cold  side  of  the  thermopile  surfaces  ( T2-T3 ),  the 
efficiency  r],  the  collector  and  radiator  area  A,  the  total  area 
of  the  thermopiles  ATE,  the  mass  of  thermoelectric  material 
Mte,  the  overall  STG  mass  ME,  the  heat  fluxes  through  the 
collector  Qs  and  the  radiator  Qr ,  the  specific  power  W/ME,  the 
power  per  unit  area  W/A2  are  reported  for  various  values  of 
ratio  K  (A/ATE). 

Taking  into  account  the  dimension  parameter,  the  results 
indicate  that  there  is  an  optimum  value  of  K  to  get  the  lowest 
area  A.  The  minimal  area  A  is  2.4  m2  matching  K  ratio  values 
between  35  and  40.  Under  these  optimal  conditions,  the 
power  per  surface  unit  reaches  167  W/m2  and  the  tempera- 
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Table  1 

Modeling  results  for  a  STG  designed  according  to  concept  I  (flat-plate),  for  different  K  =  A/ATe  values  (output  power:  400  W,  working  distance:  0.45  a.u., 
thermoelectric  legs  length:  1  cm,  insulating  thickness:  1.5  cm,  radiator  and  collector  plate  thickness:  0.1  cm,  materials  and  properties  as  defined  in  the  text) 


K 

C-t4,  T2-T3  (K) 

rj  (%) 

A(m2) 

Ate  (m2) 

A^te  (kg) 

Mf  (kg) 

Qs  (kW) 

Qr  (kW) 

W/Mf  (W/kg) 

WIA2  (W/m2) 

20 

779.3-562.8,  779.2-562.9 

3.8 

2.64 

0.132 

10.13 

36.4 

11.70 

11.30 

11.0 

151 

25 

802.9-555.3,  802.8-555.4 

4.4 

2.50 

0.100 

7.63 

32.5 

10.49 

10.09 

12.3 

160 

30 

822.1-548.6,  822.0-548.7 

4.9 

2.43 

0.081 

6.18 

30.4 

9.73 

9.33 

13.2 

165 

35 

838.1-542.5,  838.0-542.6 

5.2 

2.38 

0.068 

5.24 

29.2 

9.23 

8.83 

13.7 

167 

40 

851.6-537.0,  851.5-537.1 

5.6 

2.40 

0.060 

4.59 

28.6 

8.88 

8.48 

14.0 

167 

45 

863.1-532.0,  863.0-532.1 

5.9 

2.43 

0.054 

4.11 

28.4 

8.64 

8.24 

14.1 

165 

50 

873.0-527.4,  873.0-527.5 

6.1 

2.45 

0.049 

3.75 

28.4 

8.47 

8.07 

14.1 

163 

55 

881.7-523.2,  881.6-523.3 

6.3 

2.48 

0.045 

3.47 

28.5 

8.35 

7.95 

14.0 

160 

60 

889.3-519.4,  889.3-519.4 

6.5 

2.52 

0.042 

3.24 

28.8 

8.27 

7.87 

13.9 
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ture  difference  is  about  3 10  K.  A  relatively  high  efficiency  of 
about  5.5%  is  achieved  with  a  low  amount  of  thermoelectric 
material  (near  5  kg).  The  overall  STG  mass  includes  the 
thermoelectric  material  mass,  the  insulator  mass,  which  is 
about  8.5  kg  in  this  case,  and  the  collector  and  radiator 
masses.  The  latter  depend  of  course  on  the  AIN  plate 
thickness  of  the  collector  and  radiator.  The  influence  of 
the  thickness  (varied  between  2.5  and  0.63  mm)  on  the  key 
parameters  has  been  investigated.  It  is  found  that  the  results 
have  little  effect  except  on  the  overall  weight  of  the  STG  as 
illustrated  in  Fig.  3  for  the  best  configuration.  For  a  thickness 
of  1  mm,  the  weight  of  the  STG  is  near  29  kg  but  can  be 
reduced  down  to  23  kg  by  decreasing  the  plate  thickness  to 
0.63  mm. 

Hot  and  cold  side  temperatures  ranging  between  about 
780-890  and  520-560  K,  respectively,  are  well  in  the  range 
of  the  operating  temperature  of  skutterudite  materials.  It  has 
to  be  noted  that  the  temperature  drops  through  the  collector 
and  radiator  plates  are  less  than  0.1  K  and  can  be  therefore 
neglected. 


Fig.  3.  Dependence  of  the  overall  STG  mass  on  the  thickness  of  the 
radiator  and  collector  plate  for  various  STG  concepts  delivering  400  W  at  a 
working  distance  of  0.45  a.u.  1:  Design  I,  /  =  1  cm,  K  =  A/ATE  =  35;  2: 
design  I,  l  =  2  cm,  K  =  A/AtE  =  22;  3:  design  III,  l  —  1  cm,  K\  =55, 
Ki  =  M  /Ate  =  25  (total  mass  is  shown  without  mass  of  concentrator  (no 
collector  plate)). 


It  is  necessary  to  note  that  the  STG  efficiency  can  be 
enhanced  by  increasing  the  K  ratio  mainly  because  of  an 
enhancement  of  the  temperature  gradient  across  the  thermo¬ 
piles.  This  ratio  must,  however,  not  exceed  too  large  values 
(50  is  a  limit),  otherwise  the  transverse  heat  flows  due  to  the 
temperature  distribution  along  the  collector  and  radiator 
surfaces  will  be  considerable,  resulting  in  a  significant 
deterioration  of  the  STG  performance.  Therefore,  the  results 
for  K  values  greater  than  50  should  not  be  considered. 

Similar  calculations  have  been  carried  out  for  thermo¬ 
electric  legs  of  2  cm  in  height  (instead  of  1  cm)  coupled  to  an 
insulation  2.5  cm  thick.  The  associated  dimensions  and 
performance  resulting  from  this  study  are  summarized  in 
Table  2.  In  this  case,  the  minimal  STG  area  is  even  smaller 
than  in  the  previous  case  because  the  heat  losses  through  the 
insulating  parts  decrease  due  to  its  greater  thickness.  It  is 
equal  to  2.16  m2  at  optimal  ratio  K  =  20.  It  is  also  accom¬ 
panied  by  an  enhancement  of  the  power  per  surface  unit  up 
to  185  W/m2,  for  a  temperature  difference  of  330  K  and  with 
an  efficiency  that  reaches  5.9%.  But  consequently  the  mass 
of  required  thermoelectric  material  rises  up  to  15  kg  and  full 
STG  mass  will  also  be  significantly  greater  (about  43  kg). 
The  overall  STG  mass  is,  as  in  the  previous  case,  dependent 
on  AIN  plate  thickness  (Fig.  3). 

As  it  was  demonstrated  previously,  the  influence  of  the 
length  of  the  skutterudite  thermoelement  is  crucial  for  the 
results.  Decreasing  the  height  implicitly  means  a  decrease  of 
the  cross-section,  to  keep  their  ratio  equal  to  its  optimal 
value,  and  so  a  lesser  soldered  junction  area.  It  is  clear  that  a 
strong  reduction  of  these  surfaces  is  not  favorable  to  the 
mechanical  strength  of  the  thermopile  since  it  will  affect  the 
durability  of  the  STG  structure.  For  these  practical  reasons, 
investigation  for  leg  length  less  than  1  cm  was  not  consid¬ 
ered. 

3.2.  Concept  II:  larger  radiator  than  collector  area 

This  second  design  is  derived  from  the  first  configuration. 
It  consists  of  a  larger  radiator  than  collector  area  (Fig.  1). 
This  configuration  has  retained  our  attention  because  by 
enhancing  the  radiator  area  relative  to  the  collector  area,  one 
can  hope  to  shift  the  cold  temperature  towards  lower  values 
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Table  2 

Modeling  results  for  a  STG  designed  according  to  concept  I  (flat-plate),  for  different  K  =  A/ATE  values  (output  power:  400  W,  working  distance:  0.45  a.u., 
thermoelectric  legs  length:  2  cm,  insulating  thickness:  2.5  cm,  radiator  and  collector  plate  thickness:  0.1  cm,  materials  and  properties  as  defined  in  the  text) 


K 

Ti-T4,  T2-T3  (K) 

ri(%) 

A(m2) 

Ate  (m2) 

A^te  (kg) 

Mf  (kg) 

Qs  (kW) 

Qr  (kW) 

W/MF  (W/kg) 

WIA2  (W/m2) 

16 

837.2-542.3,  837.1-542.4 

5.2 

2.19 

0.137 

21.0 

47.5 

8.46 

8.06 

8.41 

182 

18 

850.2-537.0,  850.1-537.0 

5.6 

2.17 

0.121 

18.4 

44.9 

8.07 

7.67 

8.92 

184 

20 

861.5-532.0,  861.4-532.1 

5.9 

2.16 

0.108 

16.6 

42.9 

7.77 

7.37 

9.32 

185 

22 

871.4-527.3,  871.4-527.4 

6.1 

2.17 

0.099 

15.1 

41.6 

7.53 

7.13 

9.61 

184 

24 

880.3-523.0,  880.2-523.1 

6.3 

2.18 

0.091 

13.9 

40.7 

7.35 

6.95 

9.84 

183 

26 

888.3-518.9,  888.2-519.0 

6.7 

2.21 

0.085 

13.0 

40.0 

7.20 

6.80 

10.0 

181 

28 

895.4-515.1,  895.3-515.2 

6.7 

2.24 

0.080 

12.2 

39.6 

7.07 

6.67 

10.1 

179 

30 

901.9-511.5,  901.8-511.6 

6.9 

2.25 

0.075 

11.5 

39.3 

6.98 

6.58 

10.2 

111 

resulting  in  higher  performance.  Actually,  preliminary 
results  assuming  perfect  insulating  parts  suggested  that  this 
configuration  could  be  very  promising. 

The  calculations  have  been  made  for  various  ratios  of 
collector  to  thermopiles  areas  Kx  and  radiator  to  thermopiles 
areas  K2  assuming  a  thermoelement  length  of  1  cm.  It  has 
been  further  assumed,  for  simplicity,  that  the  surface  of  the 
insulating  parts  facing  the  Sun  has  the  same  temperature  as 
the  collector.  The  results  are  reported  in  Table  3  for  the 
particular  case  K\  =  20.  The  data  clearly  show  that  the 
increase  of  K2  leads  only  to  an  increase  of  the  STG  dimen¬ 
sions  and  mass,  and  also  to  bad  performances  (low  tem¬ 
perature  difference:  ~225  K,  low  efficiency:  ~3.8%).  These 
conclusions  remain  valid  when  increasing  K\.  The  reason  is 
that  an  increase  of  the  radiator  area  respective  to  the 
collector  results  in  an  increase  of  the  shunt  heat  transported 
through  the  insulation.  Then,  the  amount  of  heat  dissipated 
by  the  radiator  is  enhanced,  restraining  the  STG  perfor¬ 
mances.  Thus,  compared  to  the  ideal  case  (presence  of  a 
perfect  insulator)  the  real  case  is  less  attractive  and  this 
concept  may  be  discarded. 

3.3.  Concept  III:  collector  made  of  multiple  mini¬ 
concentrators 

The  last  concept  considered  is  based  on  the  use  of  solar 
concentrators.  There  are  several  types  of  concentrators  that 
could  be  used,  but  their  description  is  outside  the  scope  of 
this  paper  (see  an  example  in  Fig.  2).  The  results  of  the 
simulation  for  this  configuration  are  reported  in  Table  4  for 
the  case  K2  =  25.  It  can  be  seen  that  a  significant  decrease  of 
the  STG  area  can  be  obtained  by  increasing  the  solar 
concentration  coefficient  Kx.  For  similar  efficiency  and 


temperature  difference  the  radiator  area  is  decreased  by 
60%  as  compared  to  a  flat-plate  STG  (design  I),  resulting 
in  high  power  per  unit  surface.  It  is  clear  that  the  total  mass 
of  radiator,  insulator  and  thermopiles  will  also  be  reduced 
(note  that  MF  does  not  take  into  account  the  mass  of  the  solar 
concentrator).  This  mass  reduction,  illustrated  in  Fig.  3, 
contributes  directly  to  a  high  specific  power.  But  it  must 
be  kept  in  mind  that  the  increase  of  the  solar  concentration 
coefficient  shifts  the  hot  side  temperature  towards  the  upper 
level  and  therefore  care  has  to  be  taken  not  to  exceed  the 
operating  temperature  range  of  the  thermoelectric  materials. 
It  is  apparent  from  these  results  that  the  dimensions  and 
performance  of  this  STG  concept  using  solar  concentrators 
are  very  attractive  but  from  a  practical  point  of  view  the 
technical  development  of  such  configurations  seems  to  be 
more  complex  as  compared  to  that  of  a  flat-plate  STG. 

4.  Conclusion 

A  theoretical  simulation  of  the  performances,  dimensions 
and  weight  characteristics  of  STG  has  been  presented.  As 
active  materials,  we  used  the  potential  of  advanced  thermo¬ 
electric  materials  based  on  skutterudites.  The  analysis  of 
three  STG  designs  has  been  carried  out  for  operation  at  a 
distance  of  0.45  a.u.  from  the  Sun  and  for  400  W  require¬ 
ments.  Through  this  exploratory  work,  it  is  shown  that  an 
STG  based  on  the  use  of  skutterudite  thermoelectric  materi¬ 
als  offer  attractive  performance  features  as  primary  or 
auxiliary  power  source  for  spacecraft  in  near-Sun  missions. 
The  skutterudite  materials  are  however  still  under  develop¬ 
ment  and  need  to  be  fully  characterized  in  terms  of  perfor¬ 
mance  and  long  life  operation. 


Table  3 

Modeling  results  for  a  STG  designed  according  to  concept  II  (flat-plate),  for  K\  =  Ai/Ate  =  20  and  different  K2  =  A2/ATE  values  (output  power:  400  W, 
working  distance:  0.45  a.u.,  thermoelectric  legs  length:  1  cm,  insulating  thickness:  1.5  cm,  radiator  and  collector  plate  thickness:  0.1  cm,  materials  and 
properties  as  defined  in  the  text) 


k2 

Ti-T4,  T2-T3  (K) 

t/(%) 

A,  (m2) 

A2  (m2) 

Ate  (m2) 

A/te  (kg) 

Mf  (kg) 

Qs  (kW) 

Qr  (kW) 

W/Mf  (W/kg) 

WIA2  (W/m2) 

20 

779.3-562.8,  779.2-562.9 

3.8 

2.64 

2.65 

0.132 

10.1 

36.4 

11.7 

11.3 

11.0 

151 

30 

742.1-517.1,  742.0-517.1 

3.9 

2.60 

3.90 

0.130 

9.96 

44.7 

12.3 

11.9 

8.95 

102 

40 

713.0-486.6,  712.9-486.6 

3.8 

2.68 

5.37 

0.134 

10.3 

55.3 

13.2 

12.8 

7.23 

74 

50 

688.4-464.0,  688.3-464.1 

3.7 

2.82 

7.07 

0.141 

10.8 

68.0 

14.3 

13.9 

5.88 

56 
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Table  4 

Modeling  results  for  a  STG  designed  according  to  concept  III  (solar  concentrator),  for  K2  ^A2/AjE  =  25  and  different  Kx  values  being  here  the 
concentration  coefficients  (output  power:  400  W,  working  distance:  0.45  a.u.,  thermoelectric  legs  length:  1  cm,  insulating  thickness:  1.5  cm,  radiator  plate 
thickness:  0.1  cm,  materials  and  properties  as  defined  in  the  text) 


Ti-T4,  T2-T3  (K) 

rj(%) 

A2  (m2) 

AXe  (m2) 

37Te  (kg) 

MFa  (kg) 

Qs  (kW) 

Qr  (kW) 

W/MFa  (W/kg) 

WIA2  (W/m2) 

30 

835.9-568.9,  835.8-569.0 

4.8 

2.07 

0.083 

6.31 

20.1 

9.59 

9.19 

19.9 

194 

40 

882.7-588.1,  882.6-588.2 

5.3 

1.60 

0.064 

4.90 

15.6 

8.54 

8.14 

25.6 

250 

45 

899.8-595.2,  899.7-595.3 

5.5 

1.47 

0.059 

4.48 

14.3 

8.21 

7.81 

28.0 

273 

50 

914.1-601.1,  914.0-601.2 

5.7 

1.35 

0.054 

4.17 

13.3 

7.96 

7.56 

30.1 

294 

55 

926.2-606.1,  926.1-606.3 

5.9 

1.28 

0.051 

3.92 

12.5 

7.75 

7.35 

32.0 

312 

a  Without  mass  of  concentrator. 
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